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1. Introduction 

As has been known since Bohr, van Leeuwen and Landau edge effects play an 
important role in the physics of magnetized charged-particle systems in equilibrium. 
In particular, the diamagnetic response of these systems in the quantum regime is 
determined by electric currents flowing near the walls. The profiles of these currents, 
and of the closely related particle density, for a non-interacting magnetized electron 
gas near a hard wall parallel to the magnetic field have been investigated in detail g , 
j^l - Q. Much less is known about the profiles in an interacting magnetized electron 
gas. 

Equally important for a physical understanding of the properties of an equilibrium 
quantum system are the correlation functions. For positions in the bulk of the system 
these have been studied extensively, both for a non-interacting magnetized electron 
gas and for its interacting counterpart. For the non-interacting gas the bulk pair 
correlation function can be determined analytically both for dilute systems at high 
temperatures and for dense low-temperature systems, in which quantum degeneracy 
effects are important [Q] . For the interacting electron gas information on the behaviour 
of the bulk correlation functions is more difficult to obtain. Even for the non- 
magnetized case these functions have surprising properties. In fact, it has been 
demonstrated that the bulk correlation functions of the interacting electron gas possess 
slowly decaying tails with an algebraic dependence on the position difference || . For 
the magnetized interacting gas similar methods have been employed to prove the 
existence of analogous algebraic tails, albeit with a different exponent ||. 

The correlation functions are expected to change near a hard wall. For a 
non-magnetized free-particle system these changes are easily determined by using a 
reflection principle |l0| . The problem becomes a lot more complicated when either 
interactions between the particles or a magnetic field or both are incorporated. In a 
recent paper the interactions between the charged particles have been taken into 
account in a model system consisting of two quantum charges immersed in a classical 
plasma confined by a wall. An algebraic tail in the pair correlation function of the 
quantum particles near the wall was found. However, the exponent governing the 
algebraic tail turned out to be different from that of the bulk correlation functions 
discussed above. The result corroborates earlier findings based on linear response 
arguments . The influence of a magnetic field on the surface correlation functions 
in the adopted model system remains to be studied. 

As the influence of a wall on the correlations in magnetized quantum systems is 
not yet fully understood, it appears to be useful to try and investigate the correlation 
functions for the relatively simple case of a non-interacting magnetized electron gas 
in the presence of a hard wall. In the following we present some new results for this 
system. In particular, we will analyze the correlations for the strongly degenerate 
case of high density and low temperature, where the influence of Fermi statistics is 
important. Both strong and weak fields will be considered, so that the number of 
filled Landau levels can vary considerably. 

The paper is organized as follows. We start with two sections that serve to 
prepare the ground. In section 2 we define the relevant correlation functions for a 
system of independent particles and discuss their relation to the one-particle Green 
functions. The pair correlation function in the bulk is considered in section 3, where 
the influence of the magnetic field on the correlations is determined both analytically 
and numerically. After these preparatory sections we start considering the influence 
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of the wall in section 4. In that section we use the so-called 'path-decomposition 
expansion', which follows from a path-integral formulation, to determine the lowest- 
order corrections in the correlation functions at positions in the transition region, 
where the presence of the wall starts to be felt. An alternative way to determine these 
corrections is based on an eigenfunction expansion of the Green function, which is 
the subject of section 5. The asymptotic form of the correlation functions for large 
position differences is established in section 6, separately for directions parallel with 
and transverse to the magnetic field. In section 7 the correlation functions for positions 
close to the wall are studied, again for both directions. In the final section 8 some 
conclusions will be drawn. 



2. Correlations 

The equilibrium quantum statistical properties of a system of independent particles 
are determined by the temperature Green function 

G (r, r') = (r|e-^|r') = £ e"^" ^ B (r) <(r'). (1) 

n 

Here (3 is the inverse temperature, and ip n (r) and E n the eigenfunctions and 
eigenvalues of the one-particle Hamiltonian H, which is assumed to be independent of 
the spin of the particles. The reduced single-particle density matrix pp^r, r') of such 
a system at inverse temperature f3 and chemical potential fi is found by incorporating 
the effects of quantum degeneracy. For Fermi-Dirac particles one has 

PfiAr, r') = 2£ l + J iEn -, } Mr) <(r') (2) 

n 

where the spin degeneracy has been taking into account. The local particle density 
P/3,Ai( r ) 1S the diagonal part of (^). 

For a completely degenerate system at zero temperature the reduced single- 
particle density matrix becomes 

p M (r,r') = 2^ 0(n-E n )i, n (r)r n (r') = G^(r,r>) (3) 

n 

with 9 the step function. The diagonal part gives the local particle density /0 M (r) of 
the completely degenerate system. The ^(-dependent Green function, as defined here, 
is obtained from the temperature Green function by an inverse Laplace transform Jl3| 

i ,-C+ioo r, 

G,(r,r') = — d/3e^-Gp(ry). (4) 

27T1 Jc-ioo P 

with c> 0. 

The n-particle reduced density matrix p^(r,r') follows from its one-particle 
counterpart by a symmetrized factorization: 

n 

P ( ;l(r n y") = ]T e-n^(^.< (j -))- (5) 

Here the sum is taken over all permutations of the n position vectors, with e 71 ' the sign 
of the permutation. The structure of the n-particle reduced density matrix has been 
analyzed quite generally for a system of interacting particles by using a path-integral 
formalism Mm . The factorization property for a system of independent particles then 
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follows as a special case. The argument is not changed by incorporating an external 
magnetic field and a hard wall that confines the system. 

For the completely degenerate case a relation similar to (||) holds for (r n , r' n ) 
. In particular, the diagonal part of the two-point reduced density matrix at zero 
temperature is 



,(2) 



(rr', rr') = G p (r, r)G>', r') - G„(r, r')G^r', r). 



Often it is convenient to introduce the two-point correlation function 



g(r,r') = 



(2)/ I l\ 

Pfj,'{rr,rr') 



- 1 = -- 



(6) 



(7) 



p^r) Pll (r>) ' G^r,r)G^r',r>) 
In the following we will study this correlation function, and the influence of a magnetic 
field and a hard wall on its properties. 



3. Correlations in the bulk 



We consider a system of charged particles which move in a magnetic field directed 
along the z-axis. The interaction between the particles is neglected. To describe the 
magnetic field we adopt the Landau gauge, with vector potential A = (0, Bx, 0). The 
particles are confined to the half-space x > by a plane hard wall at x = 0. 

For positions far from the wall the temperature Green function Gp(r,r') reduces 



to the bulk Green function G 



G» (r,r>) = 



1 



B 



The latter is given by 
B 



/2¥/3 4 7 rsinh(/3 J B/2) 
x exp 



exp 



iB 

—(x + x')(y-y') 



4tanh(/3B/2) 

(z- z' f 



(r_L - rl) 5 



exp 



2/3 



(8) 



Units have been chosen such that the charge and the mass of the particles drop out, 
while h and c have been put to 1 as well. From now on we will often measure distances 
in terms of the cyclotron radius 1/yB. To that end we introduce the dimensionless 
variables £ = \f~B{x + x')/2, £ = \fB{x - x'), r\ = y/B(y - y'), and C = VB(z - z'). 
The /i-dependent Green function follows by inserting (ph in M). One finds 



G b Jr,r') = 



1 
2?ri 



c+ioo 



x exp 



df3e 



0Bv . 



B 



2 5/2 ^3/2 fjZ/2 



5o 



1 + go' 

8<Zo 



(e+v 2 ) 





[ C 2 1 


exp 


2(3B 



(9) 



with go = tanh(/3_B/4). The chemical potential is measured in terms of the energy 
difference between adjacent Landau levels, by employing the dimensionless variable 
v = p/B. Let us use now the generating- function identity |D| 



exp 



(1 



4<7o 



OO 



1 



to express the exponential function in terms of Laguerre polynomials, 
substitution of (1 — Qq) / (I + Qq) = c~P B / 2 this gives 



l - go 
l 



(10) 
After 



■ exp ( i^r/ 



X 

n=0 



c'+ic 



:— / d S e s ^-(" +1 /2)l s -3/2 
27ri ' 



e 

2s 



(11) 
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where we have set s = (3B and d = cB. The sum over n can be interpreted as a sum 
over Landau levels. 

The inverse Laplace transform in (|ll]) can be found in [ fl6| : 

i rC+lCC -i 

— / ds e st s~ 3/2 e~ a/s = -= sin (2-s/oi) Bit) (12) 
27ri J c _ ioo 

for a > and c > 0. Use of this identity in ( |ll|) results in the following expression for 
the /i-dependent Green function in the bulk: 

d3/2 / _ , „2' 



(' £±£) MV^EFllim. (13 ) 



The prime indicates that the summation is only over those values of n for which 
v — (n + 1/2) is positive, i.e. over the Landau levels that are at least partially filled. 
The bulk density follows as the diagonal part of ( |l3|) 

P\ = G>, r) = y/"-(n+l/2). (14) 

n 

The two-particle correlation function is found upon substitution of ( |l3| ) and (^) in 

©• _ _ 

The expressions (gj) and (|]J) are particularly useful for strong fields when only 
a few Landau levels are occupied. An alternative expression, which is useful for weak 
fields only, has been derived in Q. To study the limit B — > in (|l3|), we return to 
the original variables, since measuring the distances in terms of the cyclotron radius, 
or the chemical potential in terms of the energy difference between adjacent Landau 
levels, becomes meaningless for vanishing magnetic fields. The number of terms in 
the sum becomes large, as the upper limit is inversely proportional to B at fixed /i. 
Furthermore, the argument of the Laguerre polynomial gets small for fixed x — x' and 
y — y' , Hence, we can use the asymptotic form of the Laguerre polynomials |lq| 



L n (u) « e u / 2 J (V2(2n+l)u) (15) 
which is valid for u/(n + 1/2) 1 / 3 <C 1. Use of this approximation gives 
iB(x + x')(y - y')' 



B 3/2 

G>,r')«^r cxp 



J (y/2B(n + l/2)\r ± -r' ± \) 



2 

/m y2[^-B(n+\/2)] (z - z')] 
X y/B{z-z>) 

for small magnetic fields. The subscripts denote the transverse parts of the position 
vectors, which follow by projection on the xj/-planc. 

If B approaches zero, the number of Landau levels becomes very large, and their 
spacing becomes very small. Therefore it is permitted to replace the summation over 
Landau levels in (Uq) with an integral. In the limit of vanishing B we get 



G>(r,r') = \ r<ltJo(V2t|ri-ri|) 

Tf JO 



sin [V2(m -*)(*- z')] 



(17) 
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With the help of the identity JT7) 

/ dec x v+1 sin (b^a 2 - x 2 ) J v (x) 
Jo 



= J| a ^ 2 b (i + 6*)-"/3-3/4 Jv+m (a yrr^ ) 



for v = 0, we arrive at 

G»(r,r') 



7T3/2 



„/|3/2 



J 3 /a(V2?k-^l)- 



(18) 



(19) 



Note that the right-hand side is an isotropic function of the position difference, as 
should be the case for a vanishing magnetic field. We can simplify it further by using 
the explicit form for the Bessel function 



The final result is 

rib / r _/\ \ 

^ ' i n 2\ r _ r /|2 



2 sin u — u cos u 



3/2 



(20) 



cos ( y/fy \r - r'|) - — ^— sin |r - r'\) 



r — r' 



(21) 



which is identical to what one gets by starting from the temperature Green function 
for the non-magnetized system 



G h Jr,r') 



1 



(27T/3) 3 / 2 



exp 



(r 



J\2 



20 



(22) 



The 



and applying (^). The bulk density in the field-free case is = (2/i) 3 / 2 /(37r 2 ). 
two-particle correlation function in the bulk follows upon inserting (|l]) in (Q) 

In figure [l] we have plotted the bulk correlation function for B = and for B^O 
with v — 2 and v = 5. For non- vanishing magnetic field we focused on the correlation 
functions with position differences that are either parallel with or perpendicular to 
the magnetic field. For large fields, or, more precisely, for small v, the correlation 
functions for the parallel and the perpendicular directions differ considerably. For 
somewhat larger is, however, the correlation functions become fairly similar, both in 
the nodal structure, and in the amplitudes. As it turns out, these similarities are 
manifest already for v — 5, where the number of completely filled Landau levels is still 
rather low. 

Comparing (|2|) with (|l3|), we see that, by turning on the magnetic field, the 
range of the correlations in the plane perpendicular to the magnetic field becomes 
smaller, with a Gaussian instead of an algebraic decay. In contrast, the range of the 
correlations in the direction parallel to the magnetic field becomes somewhat larger. In 
fact, although the decay remains algebraic when the field is switched on, the dominant 
contribution in the tail of the correlation function becomes inversely proportional to 
the square of the distance, whereas it is inversely proportional to the fourth power of 
the distance in the field-free case. 
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Figure 1. Bulk correlation functions g(xyz, xy' z) ( ), g(xyz, xyz') ( ) 

for B ^ and (a) v = 2, (b) = 5, and g(r, r') ( — • — ) for B = 0. All curves 
start at —1 for r = r' . 



4. Path-decomposition expansion 

Introducing the wall at x = makes the temperature Green function dependent on 
the distance from the wall, i.e. the coordinate x. In the absence of a magnetic field the 
influence of the wall on the Green function is easily found from a reflection principle 
Jio| . The temperature Green function becomes 

G p (r,r') = G} 3 (r,r')-Gl(r,r") (23) 

with the bulk Green function (S) and the reflected position r" defined as 
(x",y", z") = (—x',y'z'). Likewise, the /j,-dependent Green function gets the form 

G fi (r,r') = Gl(r,r')-Gl(r,r") (24) 

with the bulk Green function 

When a magnetic field is present, the influence of the wall on the properties of 
the system is more difficult to determine. In we have seen that the temperature 
Green function can be found from a path-decomposition expansion 

oo 

G {ry) = Y,G%\ry) (25) 

n=0 

where G^\r,r') is the contribution from paths that hit the wall n times. This path- 
decomposition expansion, which was first formulated in pq |, is fully equivalent to 
the multiple-reflection expansion as introduced in pj| , and discussed for a confined 
magnetic system in Q. We have also seen in Q that for r and r' at large distances 
from the wall, the terms with small n in (|2^) are more important than those with larger 
n. In particular, the n = 1 term will give us the leading-order correction on the bulk 
quantities. The latter correspond to n = 0, so that one has G^\r,r') = G^(r,r'). If 
no field is present the expansion terminates after the term with n = 1, for all distances 
from the wall. 
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In the particular case where x — x 1 , the transverse part of the n = 1 term has 
the form 01, li 



B 2 



16tt 3 / 2 J q 



13 



& r ff)r{xy,xy') ex P 9 K ^ T {xy,xy') 



C 1 ). 



The functions f p T {xy,xy') and g'g! r (xy,xy') are given by 



f^l(xy,xy') 



{tih) 



1/2 



2 S l52(tl+t2) 1 / 2 



and 



9^l(xy,xy') = -- 



1 






4 


[(£- 





? - 2i^ 



n 



(26) 



(27) 



with ^ = tanh(rB/2), s 1 = sinh(rB/2), t 2 
sinh((/3 - r)5/2). 



h +t 2 

tanh((/3 - t)B/2) and s 2 



(28) 



4.1. Case y ^ i/ 

If we set p = q 2 /(qo 2 — q 2 ), with q = tanh[(2r — [3)B/4\ and qo = tanh(/3i?/4), we can 
write (BrJ) as 



Gy^{xy,xy') 



B 



qo 



27/2^3/2 qo 3/2 

dp 

y/p{l+p) 



exp 



2go 2 16g 

"7 1o 



tJl + (l-q 2 )p + 



2 v/l + (l-^)p 



x exp 



(l-go 2 )p ^ 1-go 2 l + (l + go 2 )p 2 ' 



(29) 



2g * 16 9o l + (l-go 2 )p 

When £ is large, only small values of (1— qo 2 )p contribute, which implies that we can set 
y/l + (1 — 5o 2 )p ~ 1- Furthermore, we can replace the last exponential function with 
exp[— (1 — qo 2 )rj 2 1 (16<7o)], at least as long as r\ is finite. If we make those substitutions 
and use 

dp 



:e- ap = e a/2 K (a/2) 

lo y/p(l+p) 

(where Kq is the modified Bessel function of the second kind), we get 



(30) 



G± <l3 {xy,xy') w - 



B 



2 7 /2 ^3/2 



t + jqo 



x exp 



4% 2 8<z ' 



go 3 / 2 
l-«o 2 



(31) 



If we are only interested in the non-degenerate case, where in general (1— qo 2 )£, 2 /qo 
is large, we can use the asymptotic expansion for Kq p5| . Multiplying the result with 
the longitudinal Green function, which is the same as in the bulk, we find the first- 
order correction to the total temperature Green function in the approximate form 

r W ( is ^ B V 1 - go 2 _ ( I 2 , -iv 

U /3 \ X V Z i X y Z )~ 27/2^3/2 ^1/2 



exp 



qo \ 2q 2 8q 

However, for the degenerate case we need the full complexity of ( ^l|) 



(32) 
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In order to obtain results for the degenerate case, we now apply (j|). In doing so 
we choose the contour of integration by setting (3 = (it + l)£/B: 



£3/2 



167T 3 ^/2 



,id7/2 



x exp 



-oo^ (it +1)3/2 
2 



t i 17 1 



l-go s 

3/2 



4g 8g ' 



1-go 2 
4go 



(33) 



In the new variables go equals tanh[(ii + l)£/4]. From £ ^> 1 one finds go « 1, which 
in turn implies 1 — rs 4exp[— £(it + l)/2]. This also means that for finite 77 we may 
replace £ + if]q /2 with £. With the help of the series representation of the modified 
Bessel function 



K (u) = J2 



"l u 
£m- 7 ~ l0g (i) 

,m— 1 



1 



2» 



and the integral relation 



dt- 



3 (it+l)a 



27TX* 



+ r(i/) 
we arrive at the asymptotic expression 



■00=) 



G${xyz,xy'z) 



r 5/2 



-f /2 p i^/2 p-^/8 V' — i 



2 2n (n!) 2 
(i/>0) 

4ri 



(34) 



(35) 



2 2 ™(n!) 2 



x 0/ - (n + 1/2) 



4[i/ - (n + 1/2)] 



n 1 



7-ln( y 



(36) 



which is valid for large vBx. Again we recognize the sum over Landau levels. In fact, 
apart from the phase e 1 ^/ 2 and the factor e - '' / 8 , the result is identical to what we 
found in flq] for r = r' . 



4.2. Case z ^ z' 

This case is considerably simpler than the y ^ y'-case, since the temperature Green 
function factorizes. The transverse part of the Green functions follows from (|3l]) by 
taking rj = 0. Using the asymptotic expansion for the Bessel function and introducing 
the longitudinal part of the Green function, one finds the total temperature Green 
function for the non-degenerate case with z ^ z' as 



G f\xyz, xy*) « ~ 27/2 J /201/2 exp (-^ - (37) 

which is the analogue of (p^). 

To calculate the /^-dependent Green function G^\xyz,xyz') for the degenerate 
case we need to keep the Bessel function in (|3l]). With s — j3B we get 

r (X)i ^ 3/2 e 1 r +i °°, 1-qp 2 ( l + go% - 2 



Stt 2 2W c _ioo * 3/2 go 3 / 2 
vA ('l^!^ex P f-^V (38) 



4go / V 2s , 

We are free to choose the contour of integration, as long as it is in the right half-plane. 
By choosing c large (i.e. c = £), we can make the same approximations for go and 
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1 — g 2 as before. In terms of the integration variable s these approximations read 
<7o ~ 1 and 1 - q a 2 w 4e~ s / 2 . Together with (pi) this yields 



G^{xyz,xyz') 



B 3/2 e -P 



2tt 2 



1 1 

X i 2im 



2 2 «(n!) 2 27rii e _ loo 

£^- ln (ir) pts 

m=l v 7 



6XP I 



(39) 



The integral in this expression is still an inverse Laplace transform. With the help of 
(O) and the analogous identity fief 



1 />C+100 

/ dse st s- 1 / 2 e- a / s 

2ni 



1 



(2Vat) 6{t) 



(40) 



for a > (and c > 0) we finally get the asymptotic expression 



^5/2 

f cm(y/2[v-(n+l/2)]C) 
\ 4[i/-(n + l/2)] 



E ^F v ^ (?1 + 1/2) 



(41) 



for large \/Bx. This result looks a bit more complicated than (p6|). In the limit £ — * 
we recover G < jj\xyz, xyz), as calculated in 

Before we discuss the expressions (H) and © in more detail, we will first show 
that the same results can be derived using a different approach. 



5. Parabolic cylinder functions 

Results equivalent to those of the previous section can be derived by using the 
explicit representation of the temperature Green function in terms of parabolic 
cylinder functions, which are the eigenfunctions of the transverse part of the system 
Hamiltonian. Because of translation invariance in the y-direction it is convenient to 
use a Fourier transform and to write the transverse part of the Hamiltonian as 

H ^ k )=-\^ + \( Bx - k ?- ( 42 ) 
For fixed k, this Hamiltonian with boundary condition ^(k,x = 0) = defines an 
eigenvalue problem. In terms of the eigenfunctions ijj n (k 1 x) and the corresponding 
eigenvalues E n (k) we define the Fourier-transformed transverse energy Green function 
1 

G±, E (k, x, x') = Mk, x)r n (k, x') S [E n (k) - E] (43) 

n 

where the eigenfunctions are normalized to dx \ip n (k, x)\ 2 = 1. In terms of this 
energy Green function the (total) temperature Green function (|l|) is 

G (r, r>) = (2^/3)~ 1 / 2 exp[-(z - z') 2 /2/3] 

x / dEe- 0E — / dke ik( - y ~ y ' ) G ±E {k,x,x'). (44) 

Jo 2lT J-oo 
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Figure 2. Correlation functions (a) g(xyz, xy' z) and (b) g(xyz, xyz') for f = 1 

( ), £ = 5 ( — ■ — ) and in the bulk (£ — > oo, ) for f = 1.5. The curves 

for f=5 and for the bulk are (almost) indistinguishable. 



m\J2(u-E)(z-z'j\ f°° ■,, , „ , s 



Performing the inverse Laplace transform as in (y) and using (|12|), we obtain 
1 

Jo 

We now choose X — X , cLS before, and switch again to the dimensionless variables £, 
£ (and a dimensionless Fourier variable n — k/y/B as well). In terms of these, the 
energy Green function is S 



Gj_^(fc, x,x) = vB 



5[E-Be n (n)]. 



(46) 



Here D\(u) is a parabolic cylinder function J15[ . Furthermore, 6„(k) is determined by 
the boundary condition 

£>e„(«)-i/ 2 (-\/2K) = 0. (47) 
If we now carry out the integration over E in (Eq) , we get 



Gfj,{xyz,xy'z') 



£3/2 



d«;e 1K?) 



sm 



(V2[i/-e„(rt)]C) 



/«„(!/) 

^ 2 „( fe )-i/2 [V2(£-«)] 



(48) 



where K n (v) is determined by e„ [k„(^)] = za 

In figure || we have plotted the resulting correlation function, which we calculated 
numerically for several values of £ and for a value of the chemical potential, which 
corresponds to a completely filled lowest Landau level. Even at moderately small 
distances from the wall the influence of the confinement on the correlation function is 
not very big. In the y-direction the correlations become a little stronger. The same is 
true for the correlations in the z-direction, at least for small £. However, the oscillating 
tail in the correlation function is suppressed by the presence of the wall. 
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Taking the limit z' — > z or £ — ► in ( |48|) is trivial. The resulting formula contains a 
phase-factor e 1Kr> , which is absent in the case y = y' as considered in We first split 
off the bulk contribution 



G^ixyz.xy' z) = - — ^ 



da e 



(K+g) ^-^+l/2 )^ ( _^ 



(49) 



The remainder G^(xy z , xy' z) — G^{xyz,xy' z) — G^xy z , xy' z) gives the correction 
due to the wall. 

In order to obtain an asymptotic expansion of G^(xyz, xy' z) for large x we split 
the integration at k' = a'£, with < a' < 1 — ^V^, and k" = a"£, with a" > \\[^- 
As in the case with r — r', the contributions from the intervals [n n {v) 1 «'] and [k" , oo) 
decay faster than exp(— £ 2 /2). With the help of the asymptotic expressions || 

e„( K )-(n + i)«^-2"e- K2 K 2 " +1 (50) 



/ d? D 2 6n(K) _ 1/2 [V2(? - k)} 
Jo 



y/nn\ 7r(n!) 2 

and 

Dl iK) - 1/2 [V2(^-K)}^rr 



2"+l g-re K 2n+l 



n , 
— ' m 



.771 — 1 



one finds 

G°(xyz,xy'z) 

with 

PnfaZ - «) = 



+ 2 2n+l (| _ K )2n ^2^+1 ^ r/^ _ x 

^2^3/2 ^ 2 2n+l 



7r(n!) 2 



v/i/-(n + l/2) 



d K e iKr > e~ K2 e"(«" K ) 2 K 2n+1 (f - «) 2n £ - k) 



4[ M -(n + l/2)] 



" 1 

E--7-ln[2^"~ K )] 



(51) 



(52) 



(53) 



(54) 



We now expand the integrand in (B3J) around k = £/2: 



dKe iK "e- K e" ( ^ K) n 2n+1 (f - nf n P n (/s,£ - k) 



we^e" 4 /2 (f/2) 4n+1 P„(C/2,C/2) / dte itn ' 2t (55) 

J — oo 

since the main contribution comes from the region around n = £/2, at least for £ 3> 1. 
Evaluating the remaining integral gives us exactly (36). 

Note that in this case, instead of G^\xyz,xy l z), we have in fact calculated 
G^(xyz,xy'z). However since the terms beyond n = 1 in the path-decomposition 
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Figure 3. Comparison between numerical data ( ) and asymptotic form 

( ) of the real part of G^(xyz, xy' z)/G'^ l (xyz, xyz) at (a)£ = 3 and (b) £ = 5 

for v = 1.5. 



expansion ( p5|) are of higher order, the asymptotic form of G c ^{xyz, xy' z) is in leading 

order identical to that of G^ {xyz,xy' z). However, in the present case we can easily 
take along more terms in ( |54| ) (see g]). The expansion around k = £/2 given here is 
only valid for rj <c £. If rj is of the same order of magnitude as £, the terms that we left 
out would not be small compared to the leading term and we would not get the right 
asymptotics. In figure |^ we compare the results of a numerical evaluation of the exact 
expression for the real part of the wall correction G^ t (xyz,xy'z)/G c ^(xyz,xyz), which 
follows from ( pl8| ) and (49), and the asymptotic form in leading order, as derived 
in (|36|). As expected, the agreement is best for small 77. The imaginary part of 
G^(xy z , xy' z) / G c ^{xy z , xy z) behaves in a similar way. 

A striking feature of the asymptotic expansion ([36|) is its Gaussian decay 
proportional to e~' ) "/ 8 . It is slower than the decay of the bulk /^-dependent Green 
function (|l3"|), which is proportional to e~ 1) / 4 . This somewhat slower decay is indeed 
corroborated by the numerical evaluation of G c ^(xyz,xy'z)/G c ^(xyz,xyz), as is shown 
in figure |[ For small 77 and large £ (i.e. for the regime where (|36| ) holds) the curves 
converge to the asymptotic value 1. 



5.2. Case z ^ z' 

On setting y = y' in ( pf ) the phase- factor e lfc,? drops out. With the help of the 
same asymptotic expressions as in the case with y 7^ y' , and the same splitting of the 
integration interval, we arrive at the asymptotic expression 

y/2B 3 / 2 ' 2 2n+1 
Gl(xyz,xy'z) « ^ ^ ^v-{n + l/2) 



7r(n!) 



f 

Jk' 



da e 



„2ri+l 



(56) 
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Figure 4. Check of the Gaussian decay, by evaluation of 

V8r]~ 1 [\og(\G'l l (xyz,xy'z)\/\G'l l (xyz,xyz)\)] 1 / 2 , for v = 1.5 and ( = 5 

( ), £ = 4 ( — ■ ■=!, £ = 3 ( ). The dotted line ( ) is the asymptotic 

value according to (p6|). 



with 

cos(,/2[^(n+l/2)R) S in(V2[^-(n+l/2)]C) 

: £; K i U — . . -. - - - - 



4[i/-(n + l/2)] v /2[^-(« + l/2)]C 

1 - 7 - In [2k(^ - «)] 1. (57) 



77) 



^ m— 1 

Again, the main contribution to the integral comes from the region around k = £/2. 
Expansion around this point allows us to evaluate the integral, and we recover ([ll|). 

In figure |^ we compare the asymptotic expression ( fill) for G c ^(xyz,xyz') with 
numerical data. As expected, the differences get smaller for increasing values of £. 

6. Correlations for large separations \r — r'\ 

The path-decomposition expansion, which we employed in section ^, is expedient for 
large x only. However, the range of validity of the representation of the previous 
section is not limited to that regime, so that we can use it for small x as well. In 
particular, it is helpful in determining the asymptotic behaviour of the correlation 
function for finite x and large distances between the points of observation, both in the 
y- and the z-direction. 

6.1. Large \y - y'\ 

Let us go back to ((48[). In order to bring the square-root singularity at k = K n (u) to 
the fore, we write 

%/2-B 3 ^ 2 x v f°° 

G il {xyz,xy'z) = ^ — / dne lKV y/K - n n (v) </>„(«) (58) 

71 n «»(") 
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Figure 5. Comparison between numerical data ( ) and asymptotic form 

( ) of G^(xyz, xyz^/G^lxyz, xyz) at (a) £ = 3 and (b) £ = 5 for u = 1.5. In 

the latter case the two curves are (almost) indistinguishable. 



with 



</>n(«) 



(59) 



12 -i/ 2 [V2(e"'-«)]' 
The function ff> n {n) is analytic on the integration interval, and proportional to 
K 2n-i/2 e -K £ or — >. qq_ Hence, the asymptotics of ( |58|) for large |t7| are determined 
by the lower boundary of the integral only. With the help of the method of stationary 
phase M we find 



K n {u) 



dne^ y/K-K n (v) n ( K ) « ^-s/V^to+^/^sg^)] r (|) 



de n (/c) 



d« 



^- 1/2 [V5((-«,M)] 



C ^ D l-m - *.("))] (60) 

where we used that the derivative of £„(k) is negative. In H we have derived the 
identity 

df^„ (re )- 1/2 (V2(f-«))' 



(61) 



This equality allows us to get rid of the normalization integral in ( pOj ) . In this way we 
arrive at the following asymptotic expression for the //-dependent Green function at 
large \n\ and finite £: 

£3/2 e i(37r/4) sgn(77) 

Gn(xyz,xy'z) ps 



2 3/2 n 5/2 



L_ V'e iK « ("to 
,3/2 e 



C/^_ 1/a [v^« n (l/)] £>"_ 1/2 [\/2(e- Kn(f))] 



dn n (v) 
dv 



-3/2 



(62) 
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Figure 6. Comparison of ri 3 / 2 ir 2 \G^(xyz, xy' z)\/(^/2B 3 / 2 ) ( j a t £ = 

v = 1.5 and its asymptotic value ( ) for large |rj|, according to (p3|) 



The asymptotic expression simplifies considerably for the special case of a 
completely filled lowest Landau level, that is for v = 3/2 and n = 0. By using the 
representation of the parabolic cylinder function in terms of confluent hypergeometric 
functions [ |l5[ we find Ko{v) — and dKQ(v)/dv = — y / 7r/2. Employing ( |6l| ) and 
inserting D\{^/2u) = \ 2u exp(— u 2 /2), we arrive at the simple asymptotic expression 

4^3/2 e i(37r/4) sgn(,,) _ 1 

Gn(xyz, xy'z) ~ ^ £ 2 e « 73^75- ( 63 ) 



for large \r)\ and v — 3/2. The asymptotic behaviour proportional to |ry| 3 ^ 2 in ( |62[ ) 
and ( p3|) is clearly induced by the presence of the wall, since the decay of the bulk 
/z-dependent Green function in the y-direction is Gaussian, at least in the presence 
of a magnetic field (see The algebraic decay is corroborated by a numerical 

evaluation of (flsf). Both results are compared in figure ^. 



6.2. Large \z — z'\ 

To determine the asymptotics for large separations in the ^-direction we set y = y' in 
(p8|). Subsequently, we need to determine the asymptotic behaviour of the imaginary 
part of the integral 



/ = 



i^2[v-e n ( K )] |C| 



n 2 



[>/§(£-«)] 



Jo~d^ n(K) _ 1/2 [x/2(f 



0] 



(64) 



for large In order to obtain the asymptotics, we split the integration interval 
at k' > 1. For k > n' we can use (50) and (|5l|). In the numerator we insert the 
asymptotic expression for the parabolic cylinder function 



^„( K )-i/2(V2(C-K))«2"/ 2+1 (-l)V™e-' t2 / 2 sinh(<~f/2), 



(65) 
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which is valid for large n and finite £. This relation can be derived by using (|0|), and 
the asymptotic relations 



D n (V2u) m (-l) n 2"/ 2 |u|"e-" 2 / 2 



and 



dD v (^/2u) 



th 



(-l) n+1 2-"/ 2 V^n!|u| 



-n-l „u 2 /2 



(66) 



(67) 



which are both valid for large negative u. In this way the contribution from n > k! to 
( p^ ) becomes 

2»-i/2 e -K 2 K 2n+1 



d K exp V2[z/- (ra + 1/2)] |C| - i- 



^nl ^Jv - (n + 1/2) 



ICI 



2 n+2 



.2" 



e"" 2 sinh 2 (/<-<f/2). 



(68) 



We now introduce a new integration variable t = e K k 2 " +1 , which implies k « \/— hrt 
for large re. The integral then gets the form 



2 n-l/2| C | 



e V2[,-(n+l/2)] Kl /"* di ( _i_ 

V^nl J Q V 0F n ! t/i/- (n+1/2) 



(69) 



sinh 2 (v^-Ei£-f 2 /2) 

x - 

-In* 

with i' = e _K re' 2n+1 <C 1. Note that there is a logarithmic singularity at t = 0. 
The asymptotic expansion of an integral of the form 
i-t' 

dte iut (-\nt)" (70) 

for t' < 1 and large u has a contribution from the lower boundary and the upper 
boundary. The contribution from the lower boundary is |2~T| 

r—k 



i E(- 1 ) r 



r=0 



E 

.fc=0 



*,r«ci> t 



(lnw)^ r 



(71) 



where r^ k '(u) is the fc-th derivative of T(u). For functions /(— hit) that can be 
expressed as a Laurent scries in v '— lnt for small i, one obtains the asymptotics of the 
integral 

rt' 

dte iut f{-lnt) (72) 

In fact, the 



(73) 



for large u by integrating the series term-by-term by means of (|7 
contribution from the lower boundary of the integral is 



u ^— ' r! 

r=0 



E 

,fe=0 



d r /(lnu) 
d(lnu) r 



Let us now return to ( |69|), which contains an integral that has the form 
of the complex conjugate of (|72|) . For this particular case one has /(—lnt) = 
sinh 2 ( v /3 Tnlf-^ 2 /2)/(-lnt) and u = 2"- 1 / 2 |C|/[^n! yjv-(n + 1/2)]. The term 
with r = in (O) dominates the contribution from the lower boundary to the 
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asymptotics, since /(In it) is much bigger than all its derivatives when u is large. 
Using this fact we get for the contribution from the lower boundary it = 0) to the 
asymptotics of (IM 



i 2n+1 c u/ 2[,-(n + i/2)l lci finh 2 (VIn^C-P/2) 



In? 



* -i V2[j-^ + l/2)] e V2[,-(n + i/2)]|ci si nh 2 (Vln^C-P/2). (74) 

The contribution from the upper boundary (t = t') to the asymptotics of ( |69|) can 
be found by use of the same theorem from [glj. It is expected to cancel against the 
contribution from the end-point k = n' of the integral 

D Lw-i /a [V2(e-/c)] 



.M Jo ^'^M-i^K'"")] 

Indeed, it is not too difficult to show that it does, by evaluating the contribution from 
the end-point at k — k' in I2 with the help of standard techniques p0[ . Finally, we 
have to check whether the end-point at k = K n (v) contributes to the asymptotics of 
I<i- The phase has a square-root singularity at that point. Using that fact, one finds 
that all terms in the asymptotic expansion that originate from the lower boundary of 
I2 are real. Hence, they drop out when taking the imaginary part of I. 
Collecting our results, we have derived the asymptotic equality 



Im/ ^ _ V2F fr + Va) ] cos (>/2[l/ _ (n + 1/2)]c) sinh2 (yfficii- e/2) (76) 

which is valid for large |£| . As a consequence we obtain the /^-dependent Green function 



25 3 / 2 sinh 2 ( VHCff- g/2) 
Gfj,{xyz,xyz ) w -5 



xJ2 y/^W ~(n + 1/2)] cos (y/2[v - (n + 1/2)] C) • (77) 



for large |C|- This asymptotic form is valid at large separations \z — z'\ and a finite 
distance 2 from the wall. A comparison with (|l^) shows that the asymptotic behaviour 
changes substantially owing to the presence of the wall. The simple algebraic tail 
(modulated by a goniometric factor) which is valid in the bulk, is replaced by a more 
subtle decay involving a logarithmic dependence on \z — z'\. 

7. Correlations near the wall 



Let us return again to (|48|). Using ( |61[) we may write it as 



Gfj,(xyz,xy z) = — —3- / d«e ^ * T |_-e„(«;) + |J 



2 - 3 r -/-m c 

x£ 2 „ (k) _ 1/2 (V2 k ) ^Lw-i/a^ - «)] (78) 

For small £ the last parabolic cylinder function can be expanded in a Taylor series. 
According to ( fl7|) the zeroth-order term of this series vanishes, so that the leading 
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term for small £ is of first order in £. The derivative of the parabolic cylinder function 
occurring in this term fulfills the Wronskian identity I15| 



D en{K) _ 1/2 {V2K)-?-D„(-V2K) , 

OK !/=e„(«)-l/2 1 L - e «W + 2 

where (|l7j) has been used again. Solving for the derivative we obtain 



^ (79) 



The right-hand side can be simplified further by introducing the integration variable 
e instead of n. In this way we arrive at the following approximate form of the in- 
dependent Green function near the wall: 



G „ ( „ V)ss ^w" dee ..,, ^°[v^K] . (81) 

„ Jn+l/2 S 

As before we consider separately the cases rj ^ and C 7^ 0. 
7.1. Case y ^ y' 

For z' — » z the Green function becomes 

o3/2 r>3/2 t2 , /-^ 

G^xyz^y'z)^ / dee iK "W" V^^. (82) 

n n Jn+l/2 

For large \rj\ the asymptotic behavior of the integral is determined by the value of the 
integrand near the upper boundary. One finds 



^3/2^/4)^) ^2 , 

_3/9i.-|3/2 2-^ e 



dv 



G»{xy Z ,xyi Zl ,~ 

The same result can be obtained from ( |62[ ) by expanding the last parabolic cylinder 
function for small f and using (f79|). 

The asymptotic decay proportional to 1 77 1 3 ^ 2 in (^3|) seems to be at variance with 
the behaviour of the correlations in the field-free case. For the latter case one finds 

4u x 2 

G M (xyz, xy'z) « y ^3 sin [ v 2 ^ (v ~ v')] • ( 84 ) 

by expanding ( p4j ) with ( ^l| ) for small x and retaining the terms dominant at large 
\y — y'\. Hence, the correlations near the wall decay faster in the free-field case than in 
the case with field. It should be noted that ( p3| ) is valid for |ry| ^ 1 or \y — y'\ S> B^ 1 / 2 . 
For B tending to the region of validity thus shifts towards 00. Furthermore, the 
number of terms in the sum becomes very large for small B (since v gets large), 
whereas the factor in front tends to 0, so that taking the limit B — ► in (B3) is not 
trivial. 

The decay of the field-free correlations can be derived in the present context by 
starting from (p^) . Let us interchange the summation and the integral and split the 
latter at e' > 1. We get 

, r u ,e' [e-V2] 

V / dee iK " (e)r ' yfV~e= \ de^V~t V e iK " (e),? 

n Jn+l/2 Jl/2 n=Q 
r u [e-1/2] 

+ / deV17~ I e iK ^ e )". (85) 

n=0 
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In the second part at the right-hand side the sum consists of many terms for each e. As 
discussed in the appendix, the values of K n (e) are located in the interval [—y/2e, V2e], 
at least for all n < [e — 1/2] — 1. As a consequence, their average separation goes to 
zero proportional to 1/ y/e. It is therefore expedient to replace the summation over n 
by an integration over a continuous variable a = K„(e)/v / 2e- The second part in ( p5| ) 
then becomes 

de^7~e / dae 1 ^" p(e,a) (86) 



with p(e, a) da the number of values of n n (e) in the interval [v2eer, v 2e(<r + da)]. For 
large \r)\ the dominant contribution in the integral over a comes from the endpoints. 
Since the density p(e,a) at the endpoints a = ±1 is 2 3 / 2 e(l =F a) 1 / 2 /it (see the 
appendix), the expression ( |S6| ) is for large \rj\ 

o3/4 rv 

/ deV^ e V4 cos (V&M-37T/4). (87) 

The integral over e can likewise be evaluated for large as once again only the 
endpoints of the integral contribute. The upper boundary gives 

-V2^ sm(V2vr)) (88) 

If 1 

while the contribution from the lower boundary in (|87j) will not be needed. 

As a final step we have to consider the first part in ( |S5| ) . Here the summation can 
be replaced by an integration for e ^> 1 only. Of course, the contribution of the upper 
boundary to the asymptotic expression has to cancel that of the lower boundary in 
(|S7|), since the final result should not depend on the choice of e'. The lower boundary 
of the integral over e in the first part of (|||) does not contribute either. In fact, only 
the term with n — survives for e near 1/2 and the factor Ko(e) goes to infinity for 
e — ► 1/2, so that for large \r/\ the contribution from the lower boundary damps out by 
interference. 

Collecting the results, we have found that for B — > and large \r/\ the asymptotic 
behaviour of the /^-dependent Green function near the wall is given by d82| ) with ( p8[ ) 
inserted: 

Gu{xyz,xy'z) w r-^ — sin (v2f rn. (89) 

it A rf 

After restoring the field-independent variables this result coincides with (|84|) . 
7.2. Case z ^ z' 

In this case the dependence on K n (e) drops out from (pl|): 



G^xyz,xyz ) » — > / de L - (90) 

W „ in+l/2 t, 

The integral is elementary. As a result we get for the Green function at small x 
o R3/2 p , , 

G»(xyz,xyz')* ^~r~ E { VW 1 ~(n+ 1/2)] C cos [y/2[v - (n + 1/2)] C] 

^ n 

-sin[V2[i/-(n + l/2)]C]}. (91) 
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By taking the limit ( -> we obtain a simple form for the particle density near the 
wall: 

o5/2 R 3/2t2 , 

^ = G^(r,r)« W-(n+l/2)} 3/2 (92) 

n 

which should be compared to the bulk density jljj). 

For large separations the first term between the curly brackets in ( |9l| ) is dominant, 
so that the /i-dependent Green function for large separations and small distance from 
the wall becomes 

OZ33/2 ?2 , 

G^xyz, xyz') « ^2^E V2[f - (n + l/2)]cos [v^FF+WC] • (93) 



The same expression can also be found from (77), by putting £ -C 1 and expanding 
the hyperbolic function for small values of its argument. 

For small B it is convenient to return to the original variables x, z — z' , fi, and 
to introduce the new integration variable u = eB in (|9^). Replacing in addition the 
summation by an integration we get 



«K f * r d „ si °'v«<— a . (94 , 



7T 

Performing the integrations we obtain the expression 

,.2 



G fl {xyz ) xyz / ) » ~ ^ 2(z _ z , )3 ( 2sin [v 2 M ( 2 - 2 ')] + cos [v 2 M (z - z' 



^5 sin [v^I («-«')]} 



(95) 



/i(z — z 

for the /x-dependent Green function at positions near the wall in the field-free case. 
It can easily be checked that the same expression is found by expanding the general 
form (|24|) for small x. 

At large separations \z — z'\ the dominant term in ( p5| ) is the first one. It agrees 
with (p4), when y and z are interchanged. 



8. Discussion and conclusion 



Our results show that the presence of a magnetic field and of a wall leads to remarkable 
changes in the pair correlation function of a completely degenerate non-interacting 
electron gas. In the bulk the correlation function follows from (0), with the in- 
dependent Green functions ( |l3|) and (H|) for the cases with and without field. Whereas 
the correlation function is isotropic and decays algebraically (cx r~ 4 ) in the field-free 
case, it becomes anisotropic with a Gaussian dependence in the transverse direction 
and an algebraic one (cx r~ 2 ) in the parallel direction. 

In the neighborhood of a plane hard wall the pair correlation function becomes 
anisotropic even for the field-free case. Its form, which follows directly from a reflection 
principle, is given by (0) with (pi|). For two positions far apart, but at equal distance 
from the wall, the ensuing correlation function decays algebraically, as in the bulk. 
However, the decay is proportional to r~ 6 , as follows from (|84|). 

If in addition to the wall a magnetic field is present, the pair correlation function 
in the vicinity of the wall becomes anisotropic for two reasons, as both the field and 
the wall break the symmetry. The general expression for the correlation function, for 
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arbitrary field strengths and arbitrary distances from the wall, follows by substitution 
of pq ) in (Q). Relatively far from the wall the corrections to the bulk correlation 
function are still small. The first-order correction terms were given in 

(§|) and @, 

and checked numerically in figures |^ and ^. Near the wall the correlation function - 
and in particular its tail for large separations - is modified considerably as compared 
to its form in the bulk. In fact, the qualitative difference between the behaviour in 
the directions parallel and transverse to the field, which is a prominent feature of the 
bulk correlation function, is lost in the vicinity of the wall. In both directions the tails 
become algebraic, albeit with a different exponent. This is seen by inspecting ( |62|) 
(or (|8^)) for the transverse direction and ((7?]) (or (|93|)) for the parallel direction. In 
the former case the decay is proportional to r -3 , whereas in the latter it is oc r~ 4 . 
Qualitatively, the change in the decay of the transverse correlation function from 
Gaussian in the bulk to algebraic near the wall can be understood in a semiclassical 
picture. In the bulk the cyclotron motion of the particles leads to a strong localization 
of the correlations, which is associated with a Gaussian decay. On approaching the 
wall the so-called 'skipping orbits' along the wall become important. They lead to 
a derealization effect in the particle motion, which manifests itself as an increase in 
the range of the correlations. In this way, the cross-over to an algebraic decay of the 
correlation function finds an explanation. 

In conclusion, we have shown that for a non-interacting electron gas the influence 
of a wall on the correlations is quite considerable, especially in the presence of a 
magnetic field. It would be interesting to determine the edge effects in the correlations 
for an interacting electron gas in a magnetic field. 



Appendix A. Zeros of the parabolic cylinder functions for large argument 
and index 



The zeros of D\{z) satisfy the inequality \z\ < 2^J\ + 1/2, with the possible 
exception of a single negative zero fl22fl . For A ^> 1 the zeros in the neighborhood of 
\z\ — 2^J A + 1/2 can be determined by approximating the parabolic cylinder functions 
by Airy functions. 

For positive z and a — z/[2y / A + 1/2] slightly smaller than 1, one has |^2| 

D x (z) « 2V2+V3 r(|A + |)(A + i) 1/6 (73^2 ) V4 Ai ([4( A + W^) (A.l) 
with r given by 

2/3 



{f[iarccos(a)-iayi^]} « -2~^{1 - a). (A.2) 



For x 3> 1 the Airy function may be approximated as |2^ 

1 



Ai(-z) 



/4 sin(§z 3 / 2 + f 



(A.3) 



Using this approximation in (A.l), one finds the zeros of the parabolic cylinder function 
from the zeros of the sine function. In terms of the variable <r, the zeros near a = 1 
are found as 

"3tt(to- 1/4)1 2/3 



1 



1 



2(A + l/2) 



(A.4) 
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with positive integer m. Since A is large, the zeros are closely spaced. In view of the 
main text we introduce e = A + 1/2 instead of A. Writing the number of zeros between 
a and a + da as p(e, a) da we get 

p(e,a) = — e(l-a) x / 2 (A.5) 

7T 

as the density of zeros near a = 1. 

Similarly for negative z and a slightly bigger than —1 the parabolic cylinder 
function can be written as p2[ 

D x (z) » 2V2+V3 r(|A + i)(A + i) 1 / 6 V4 

x {cos(^A) Ai ([4(A + i)] 2/3 r) - sin(TrA) Bi ([4(A + ±)] 2/3 r) } (A.6) 



with r w -2 _1 / 3 (1 + cr). With the use of (|A.3j ) and the analogous relation 

1 



Bi(-x) 



,1/4 



(f ^» + f ) 



(valid for large x), one obtains the zeros of D\ (z) near a = — 1 as 



-1 



1 



3vr(m + A - 1/4) 



2/3 



2(A + l/2) 

for integer m with m > —A + 1/4. The density of the zeros is found to be 

2 3/2 
7T 

for a near — 1. 



p(e,a) 



■e(l + a) 1 / 2 



(A.7) 



(A.. 



(A.9) 
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